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Abstract 

We consider the interaction of a nonlinear Schrodinger soliton with a spatially 
localized (point) defect in the medium through which it travels. Using numerical 
simulations, we find parameter regimes under which the soliton may be reflected, 
transmitted, or captured by the defect. We propose a mechanism of resonant energy 
transfer to a nonlinear standing wave mode supported by the defect. Extending 
Forinash et. al. [1], we then derive a finite-dimensional model for the interaction of 
the soliton with the defect via a collective coordinates method. The resulting system 
is a three degree-of-freedom Hamiltonian with an additional conserved quantity. We 
study this system both numerically and using the tools of dynamical systems theory, 
and find that it exhibits a variety of interesting behaviors, largely determined by 
the structures of stable and unstable manifolds of special classes of periodic orbits. 
We use this geometrical understanding to interpret the simulations of the finite- 
dimensional model, compare them with the nonlinear Schrodinger simulations, and 
comment on differences due to the finite-dimensional ansatz. 
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1 Introduction 



In a previous study, involving the first and last authors [2], (Bragg) resonant 
nonlinear propagation of light through optical waveguides with a periodically 
varying refractive index profile and localized defects was investigated. In that 
work, an approach to the design of spatial defects in a periodic structure 
for the purpose of trapping and localizing light pulses was suggested and ex- 
plored. The technique involves resonant transfer of energy from traveling waves 
(gap solitons) to nonlinear standing wave modes localized at the defect. The 
cubic nonlinear Schrodingcr (NLS) equation with localized potentials, which 
we study in the present paper, provides a simpler model exhibiting similar 
phenomena that is more amenable to analysis In particular, in numerical sim- 
ulations of NLS solitons incident on a single delta-well (point) defect in a 
one-dimensional continuum, we find a variety of behaviors depending on the 
parameters describing the soliton. 

Several studies have examined the propagation of nonlinear waves through 
variable or random media. The approach taken in [3,4] is to view such a 
medium as sequence of individual weak scatterers, each modeled by a repul- 
sive delta function potential barrier. The interaction of a soliton with an indi- 
vidual scatterer is formulated as a mapping of internal soliton parameters; a 
soliton which interacts weakly with a scatterer adjusts its internal parameters 
slightly due to radiative loss of energy. The interaction with the full medium 
is approximated by repeated composition of this simple mapping. 

The problem we address is different in a number of respects. Wc consider 
strong interactions of a soliton incident on a single scatterer or defect poten- 
tial. Furthermore, our potential is taken to be an attractive delta function 
potential well, which has a single localized eigenstate (defect mode). There- 
fore, an incident soliton can be expected to break up into a soliton with ad- 
justed parameters, due to energy transfer to the localized defect mode, and to 
outgoing radiation. Components of the soliton's energy may be refiected by 
the potential, transmitted through the potential, or captured by its intrinsic 
modes. These strong nonlinear scattering interactions exhibit a great deal of 
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complexity, which we explore first by direct numerical simulation and then via 
finite dimensional models. 

In particular, we first conduct a series of numerical experiments on the partial 
differential equation (PDE), in which a variety of of phenomena are observed, 
which we may partly explain by a resonant transfer of energy to standing wave 
modes localized at the defect. Second, we derive a finite-dimensional system 
of ordinary differential equations (ODEs) that models the interaction of soli- 
tons with nonlinear standing wave 'defect' modes supported by the potential. 
This part of the analysis is similar in spirit to our earlier study of a finite 
dimensional reduction of the simpler case of kinks interacting with a trapped 
mode in the sine-Gordon equation with a point defect [5]. After reviewing the 
basic PDE model in Section sec:pde and 3, and describing the results of direct 
numerical simulations in Sections 4, we outline the (formal) finite-dimensional 
reduction procedure in Section 5. We then describe in Section 6 a representa- 
tive set of three numerical experiments, in the same parameter ranges as the 
PDE studies, that reveal the kinds of soliton transmission, reffection and tran- 
sient capture behaviors that the ODEs exhibit. Section 7 is devoted to analysis 
of the ODEs. We describe invariant subspaces and special sets of orbits, fo- 
cusing on the stable and unstable manifolds of certain periodic orbits. These 
are shown to partially 'organize' the global dynamics, in particular provid- 
ing separatrices between transmitted and reflected soliton orbits. In the flnal 
Section 8, we make comparisons between the PDE and ODE dynamics and 
summarize. Detailed analytical calculations, and some background material, 
are relegated to an Appendix. 

In addition to our specific study of NLS soliton-defect mode interactions, we 
believe that the detailed comparison of PDE and ODE solutions in this paper 
has general implications for similar 'collective coordinate' finite-dimensional 
representations commonly used to study dynamical interactions of continuous 
fields. 



2 The nonlinear Schrodinger (NLS) equation with a point defect 

We consider a nonlinear Schrodinger (NLS) equation with a spatially localized 
'attractive' impurity (a defect potential well) at the origin: 



In the absence of a defect (7 = 0), this system supports a two-parameter 
family of solitons of the form 
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Fig. 2.1. The frequency and amphtude of a sohton in the absence of a defect (sohd), 
and of the defect mode (dashed). 

where u = — 1(?7^ — v"^) is the temporal frequency. Sohtons are nonhnear 
bound states which play a fundamental role in the unperturbed (7 = 0) NLS 
equation, and we are especially interested in their behavior in the perturbed 
system. For 7 > 0, (2.2) does not solve (2.1), but far from the defect, the 
soliton propagates essentially without distortion and at constant speed v due 
to the exponentially small overlap of soliton and potential. 

Equation (2.1) also supports exac^ nonlinear bound state or defect mode solu- 
tions of the form 

UBef{x, t) = ae*"'*/^ sech {a\x\ + tanh~^ -), (2.3) 

for all a > 7. These solutions are constructed from a stationary {v — 0) sohton 
on each side of the defect pasted together at x = to satisfy the conditions of 
continuity at a; = and the appropriate jump condition in the first derivative 
at a; = 0, m(0+) — m(0~) = —27^(0). For both bound state families, usoi and 
MDef, the frequency of oscillation depends nonlinearly on its amplitude. 

The chief concern of this paper is to understand strong interactions between 
solitons and the delta well defect. If the defect strength 7 is small of the soliton 
velocity is large, then the interaction is weak: a small amount of energy is lost 
to radiation, and the soliton continues past the defect with minor changes to 
the parameters that define it. In the case of weak interactions, estimates for 
the energy loss of the soliton can be obtained by the first order perturbation 
theory (the Born approximation) [4,3]. When 7 is sufficiently large and v is 
small, stronger interactions can take place, the character of which may be 
understood in terms of a nonlinear resonance that in some cases takes place 
between the soliton and the defect mode. 

Solitons with v = have ||Msoi||i2 = '^V and frequency —rf/2, whereas non- 
linear defect modes have HiiDefllia — 2(a — 7) and frequency — a^/2. In Fig- 



ure 2.1 we plot the squared norm of these two types of mode as functions 
of frequency. In the following section we discuss the bifurcation diagram of 
Figure 2.1 and its implications, and review the well-poscdncss theory of the 
initial value problem and the dynamic stability theory for u^oi and woef- Due 
to the relation between the square of the norm with the electromagnetic 
energy in the context of optics, wc shall refer to H^soiHiz and ||MDef|||2 as the 
energy of the soliton and defect modes, and more generally to the square of 
the norm of u over a region of space as the energy contained in that region. 

Note that no defect modes exist in the frequency range lo G (—7^/2, 0]. This 
observation is crTicial in predicting which solitons will be trapped and which 
reflected by the potential. We find, roughly, that sufficiently slow solitons with 
77 > 7 are trapped upon encountering the defect, while slow solitons with r] < j 
are reflected by the defect. This suggests that trapping occurs via resonant 
energy transfer from the soliton to the defect mode. If the incoming sohton 
has frequency less than —7^/2, it can and may excite a nonlinear defect mode 
and transfer its energy to that mode. Otherwise, as in [6], wc find that the 
defect behaves as a scatterer, splitting the incoming wave into transmitted, 
captured, and reflected parts. This behavior differs sharply from one seen in 
nonlinear optics: in the nonlinear coupled mode equations describing the in- 
teraction of gap solitons with nonlinear defect modes [2], it was found that 
in the absence of resonance, pulses were either coherently reflected or trans- 
mitted after interacting with the defect, with little energy captured or lost to 
radiation. 



3 Overview of well-posedness and stability 



Structural properties of NLS: The nonlinear Schrodinger equation (2.1) is a 
Hamiltonian system, which can be written in the form: 

iut = —n[u,u*], (3.1) 
where H[u,u*] denotes the Hamiltonian: 

■ 75(x)|ii|^ ^ dx 

dx- 7|'u(0)|2 . (3.2) 
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Invariance with respect to time-translations implies that 7i[u,u*] is conserved 
by the flow generated by (3.1). Additionally, invariance under the transforma- 
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tion u I— > e'^u, ^ & implies that 

Af[u,u*] = J\u\'^dx= \\u\\l (3.3) 

is a conserved integral. 

For the spatially translation invariant case, 7 = 0, NLS has the Galilean 
invariance: 

u{x,t) ^ u{x-vt,t) e'^''''-'^''''\ V e R\ (3.4) 

Well-posedness theory: The functional TC[-] and Af[-] are well defined on 
H^{R^), the space of functions /, for which / and d^f are square integrable. 
It is therefore natural to construct the flow for initial data of class H^. In fact, 
it can be shown that, for initial conditions uo = u{x,t = 0) e H^{R^), there 
exists a unique global solution of NLS, u e C°(i?^; H^{R)), in the sense of the 
equivalent integral equation: 

u(t) = U(t)uo +i tuit-s) \uis)\^uis) ds, (3.5a) 
Jo 

U{t) = exp {-iHt) , H = ~dl - -f5{x). (3.5b) 

The spectral decomposition of B. is known exphcitly [7] and can be used to 
construct C/(t) explicitly. 

To show the existence of a solution to (3.5a) in if ^, we must show the existence 
of a C^{R^; H^{R)) fixed point of the mapping u{x,t) 1— > J[u\{x,t), given by 
the right hand side of (3.5a). 

We now outline the key ingredients of the proof. To bound J[u] and its first 
derivative in L^, we introduce the operator A — I + PcH, where Pc denotes the 
projection onto the continuous spectral part of H. Note that ^ is a nonnegative 
operator, since the continuous spectrum of H is the nonnegative real half-line. 
Moreover, we expect ||^2/|| ~ WIWh^ = II (-^ ~ •9^)^/11^2. In fact, we shall use 
that the following operators are bounded from to L^: 

This follows from the boundcdncss of the wave operators on [8]. Therefore, 
we have an equivalence of norms 

Ci 11/11^,1 < ||^^/|U2 < C2 11/11^1. (3.6) 

Our formulation (3.5a) and introduction of A is related to the nice prop- 
erty that A, and hence also functions of A, commute with the propagator 
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exp{—iHt). We shall also use the Sobolev inequality: 

\f{x)\'< C||/|U.|I^./IU^, (3.7) 

and the Leibniz rule [9]: 

\\ii-dl)Hf9)\\ < c{\\fh^ W-dDhh^ + W - d'Jfh^ Ml^). 

(3.8) 

Since U{t) is unitary in L^, we have 
\\Ah[u]\k. 

< WA^UoWl^ +7 / \\A^u{s)\\{s)\\l2 ds 
Jo 

= WA^uoh^ + jJ^\\{AHl-dl)-'^)-{{I-dl)^u{s)\Ms))\\L^ds. 

(3.9) 

By (3.6)), (3.7), and (3.8), 

||JM(i)||^i < C||^JM|U2 < CiWuoWm + C2T sup \\u(s)\\%,. (3.10) 

se[o,T] 

Now assume that u is such that sup^g^ j^] ||-u(s) < 2Ci. Then, by (3.10), 
by choosing T < Ti sufficiently small. sup5g[o^r] || J[ti](s)||iji < 2Ci||tio||//i. 

\\J[u]it)\\m < Ci\\uo\\h^ + C2T (2Ci||«o|ki)'. 

It follows that for < T < Ti sufficiently small, the transformation J[-] maps 
a ball C°([0, T]; if^(i?)) into itself. A similar calculation shows that 

\\J[u]{t) - J[v]{t)\\m < KT {2C4u4H^f sup \\u{s) - v{s)\\h^ , (3.11) 

se[o,T] 

and therefore for < T < T2 < Ti, the transformation J[-] is a contraction 
on this ball. Therefore, J[-] has a unique fixed point in C°([0, T]; iJ^(i?)) for 
T sufficiently small and local existence in time of the fiow follows. Global 
existence in time follows from the d priori bound on the norm of the 
solution implied by the time-invariance of norm and Hamiltonian. 

Nonlinear hound states: Bound states are an important class of solutions hav- 
ing the form: 

Uh{x,t) = e-'^^ip{x;X), cpeL'^. (3.12) 

For the linear Schrodinger equation, the functions (p are eigenstates of a 
Schrodinger operator: —^d^ — 'jS{x) and satisfy the equation 

-Ti^xx - 'yS{x)(fi = X(f. (3.13) 
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Bound states are known to play a fundamental role in the general dynam- 
ics of the linear Schrodinger equation. This is a consequence of the spectral 
decomposition of linear self-adjoint operators. 



For NLS, such bound states satisfy the equation: 

~<fxx - (l^ix) + \<p\^)<P = A(^, (3.14) 

and have the general character of 'nonhnear eigenstates', although there is 
no rigorous decomposition theory of solutions into such states, except in the 
completely integrable case 7 = [10]. For this translation-invariant case there 
is family of solitary traveling wave solutions (2.2). These are Galilean boosts 
of the basic solitary standing wave (2.2) with v = 0; see (3.4). For 7 > 0, the 
equation is no longer translation-invariant and we have the defect or 'pinned' 
states of (2.3). These two famihes of nonhnear bound states are plotted in bi- 
furcation diagram of Figure 2.1. We note that for 7 = the family of solitons 
bifurcates from the zero state at zero frequency, the cndpoint of the continu- 
ous spectrum of the linearized operator —d^ about the zero state. For 7 > 0, 
the family of defect states bifurcates from the zero state in the direction of 
the eigenfunction of the linearized operator, —d^ — "y5{x), and at the corre- 
sponding eigenfrequency A = —7^/2; see [11] for a general discussion. 

Stability of nonlinear bound states: An alternative characterization of the non- 
linear bound states ?/soi and uocf is variational. The advantage of the varia- 
tional characterization is that it can be used to establish nonlinear stability 
of the 'ground state; see [12,11]. 

Theorem 1 (I) The families of nonlinear bound state profiles rj ^ (fisoii^', rj) 
for the case 7 = and a 1-^ ipoefix] a) for the case 7 > can he characterized 
variationally as minimizers of the Hamiltonian, 7i, subject to fixed L? norm, 
M: 

mm ny\, My\ = p (3.15) 

Thus Msoi and WDef are called ground states in their respective cases. Their 
associated frequencies, A(p) = —rf/'i and A(p) = — a^/2 arise as Lagrange 
multipliers for the constrained variational problem (3.15). As p \{p) 0, 
respectively, X{p) —7^/2. 

(11) Ground states are nonlinearly orbitally Lyapunov stable, i.e. if the 
initial data are close to a soliton (modulo the NLS symmetries), then the 
solution remains close to a soliton in this sense for all t e (—00, 00). 

For results on asymptotic stability of nonlinear ground states, see [13-17]. 
High energy defect modes and solitons: Solitons of the translation-invariant 
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NLS equation {usoi of (2.2) for 7 = 0) may be related to the high energy (L^ 
norm) nonhnear bound states, WDcfj in an appropriate limit. The following 
summary is based on the variational principle of Theorem 1, and does not 
require the exphcit formula (2.3) for the defect mode. This argument can also 
be applied to more general linear potentials y(a;)|wp in the Hamiltonian in 
place of 5(x)|iip, and to more general nonlinearities. 

Consider the variational problem (3.15), in which we make explicit the de- 
pendence of the Hamiltonian Ti on the defect strength 7 by writing 7i{u; 7) 
in place of the notation of (3.2). Define Tp[M](a;) = pu{px), for p > 0. Then, 
-^T-i{Tp[u]; 7) = T-C{u; p. If I{p; 7) denotes the minimum in (3.15), we there- 
fore have: 

;^^(p;7)=:^(i;^)- 

Note that, as p ^ 00, X(l; ^) formally approaches X(1;0), the constrained 
minimum of the cubic NLS Hamiltonian 7i(-u; 0), and that the extremizer is the 
classical one-soliton. It can be shown that if u-De{{x; a) denotes the nonlinear 
defect mode and we define Ud hy Ti/a[uDef]{x;a) = UD{x;a), then Ud^x^o) 
converges strongly to the classical one-sohton of norm one. It follows that, 
for large a, Ujieii^, 0) looks more and more like a Ta-scaled soliton (a solitary 
standing wave) of the translation-invariant NLS. 



4 Direct Numerical Simulations of the PDE 

In this section we discuss simulations of the initial value problem for the NLS 
equation (2.1). All numerical experiments in this section were performed using 
a modification of a finite difference approximation due to Fei, Perez-Garcia 
and Vazquez [18], which conserves a discrete norm and, in the absence 
of a defect, a discrete analog of the Hamiltonian. The method is accurate to 
second order in both space and time, and is implicit only in its linear terms. 
Therefore, it requires a linear equation to be solved at each step. The Dirac 
delta function is approximated either as a single point discontinuity, or by a 
smoother function with very small support, with similar results in either case. 

In the numerical experiments, a soliton is initialized far from the defect with 
prescribed velocity v and amplitude 7] and is allowed to propagate toward 
the defect location. A wide variety of behaviors is seen as parameters are 
varied. For simplicity (and motivated by a scaling argument given in section 5) 
we limit our study to defects with strength 7 = 1. Therefore, the branch 
of nonlinear defect modes bifurcates at the frequency —1/2. The nonlinear 
resonance summarized in Figure 2.1 is useful in understanding the various 
behaviors that are possible in this interaction. In the figures that follow, \u\ 
is plotted although the system's conservation might suggest plotting \u\'^. 
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Fig. 4.1. Input versus output velocities of a soliton with rj = A. 




Fig. 4.2. The soliton amplitude \u\ after interaction with the soliton, with large initial 
soliton amplitude ij = A. On the left, a slower soliton is captured. On the right, a 
faster soliton is transmitted, leaving behind a small defect mode. The transition 
from capture to transmission is abrupt, occurring at a critical velocity Vc ~ 1.78, as 
seen in Figure 4.1. 

This is done to render visible the radiation that is shed during the interaction. 

In the first set of runs, rj is set to 4, and v is varied between 1 and 2. In this 
case, the solution remains mainly soliton-like, with very little loss of energy 
to radiative modes. There exists a critical velocity, above which the soliton is 
transmitted past the defect, although with diminished speed, and below which 
the soliton transfers its energy to a nonlinear defect mode. Figure 4.1 shows 
the input versus output velocities, indicating the critical velocity v^. ^ 1.78. 
Figure 4.2 shows the evolution of |M(a;, t)| up to a time somewhat after the 
interaction has taken place. In the left hand figure, the initial velocity was 
about V — 1.65, and after some time the solution is a defect mode centered 
at a; = with a small amount of radiation. In the right hand figure, with 
initial velocity greater than Vc, the soliton largely survives, although some of 
the incoming soliton's energy is captured by the defect and, as time proceeds, 
eventually takes on the form of a small amplitude defect mode. 
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Fig. 4.3. A soliton with initial amplitude ij = 2 has a large amount of its energy 
captured by the defect. At left, a slower soliton with v = 1.5 has a substantial portion 
of its energy captured. At right, as the incoming soliton's velocity is increased to 
-y = 1.75 , less energy is captured and a larger soliton gets through. 
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Fig. 4.4. For 77 = 2, as the input velocity is increased, the total amount of energy 

captured by the defect mode decreases. 



When 7] is reduced to 2, the behavior changes. The simulations show a more 
complicated nonlinear scattering process. The pulse sphts into three parts: 
reflected, captured, and transmitted. In all cases, a significant defect mode 
is created. The faster the soliton's initial speed, the smaller the defect mode 
remaining in a neighborhood of the origin and the larger the transmitted por- 
tion. An example is shown in Figure 4.3 and the phenomenon is summarized 
in Figure 4.4, which shows how the fraction of energy deposited into the defect 
mode decays monotonically as the input velocity increases. 

When 1] is further decreased to 0.5, the behavior alters yet again. In this case, 
only a small amount of energy is captured by the defect, while large amounts 
are reflected and transmitted; see Figure 4.5. When the incoming velocity is 
sufficiently small, the soliton appears to be completely reflected and when 
the incoming velocity is sufficiently large, the soliton is nearly completely 
transmitted. For intermediate velocities, the pulse is split into a transmitted 
and a reflected wave. This nonlinear scattering phenomenon is studied by 
Cao and Malomed [6], who derive approximate reflection and transmission 
coefficients for the interaction in the case of small 77. 



11 



Fig. 4.5. When rj = 0.5, most of the wave is reflected with v = 0.25 (left) or 
transmitted when v = 4 (right). When t; = 1, the wave is spht nearly in half into 
reflected and transmitted portions. 

4-1 Discussion of Simulations 

We wish to interpret the scattering of sohtons by defects in terms of the 
amplitude- frequency curves of Figure 2.1. In the previous section, as the soliton 
amplitude parameter rj is decreased, two effects make the transfer of energy 
from the soliton to the defect mode less efficient and enhance the scattering. 
At first, for 77 = 2, the amount of energy lost in the interaction is increased, 
compared to the experiments with r) — A. Finally, for 77 < 7, there no longer 
exists a nonlinear defect mode that resonates with the incoming soliton, and 
therefore almost no energy is captured. 

It was suggested in [2] , in the context of the nonlinear coupled mode equations, 
that for sufficiently slow incident solitons, a simple resonant energy transfer 
mechanism should hold. In particular, a good approximate predictor of the 
distribution of trapped energy would be given by the vertical projection of 
the point corresponding to the incident soliton onto the corresponding point 
on the defect mode curve with the same frequency. It turns out that this 
approximation is valid in certain cases, but extensive further simulations have 
shown the general situation to be more complicated. In this subsection we 
explore this issue by means of the following auxiliary numerical experiment. 

We initialized a family of solitons with zero velocity and varying amplitudes 77, 
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Fig. 4.6. A soliton initialized over a defect quickly decomposes into a defect mode 
plus radiation. 
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Fig. 4.7. Arrows connect soliton initial conditions to defect mode steady periodic 
solutions. 

centered directly over the defect, and let them run forward until they formed 
standing wave states (recall that the defect standing wave (2.3) is an exact 
solution of (2.1)). The solutions rapidly evolved into a combination of nonlin- 
ear defect modes and radiation. If the above picture applied, then the defect 
modes thus produced should have almost the same frequencies as the initial 
solitons. An example is shown in Figure 4.6, which was initialized with rj — 1 
and 7=1. Figure 4.7 reproduces the two curves of Figure 2.1 with arrows 
connecting each soliton's initial conditions to the periodic solution of the cor- 
responding defect mode following the interaction. The arrows are far from 
vertical, and show a consistent downward frequency shift. Evidently, in this 
experiment the above 'direct' resonant exchange mechanism does not apply. 

However, this finding does not contradict the role of resonant energy transfer 
in trapping an incoming soliton. In the above experiment, in which the 'inci- 
dent' soliton has velocity v = 0, the mechanism differs from that involving a 
moving soliton. Were the nonlinear term absent, the results of this auxiliary 
experiment would be well-understood. Spectral theory dictates that the initial 
condition decomposes into a bound state and radiative modes, and there is no 
role in this interaction for the soliton's internal frequency. In the presence of 




13 



nonlinearity, the behavior is essentially the same, although we can no longer 
compute the captured solution as the projection of the initial conditions onto 
the bound state. The frequency is shifted because the solution is in some 
sense finding the nonlinear projection of the initial condition onto the family 
of nonlinear defect modes (2.3), which has only a shghtly smaller norm. 
As a is increased, the energy (L^ norm) is increased and, as discussed at the 
end of section 3, the ground state defect mode approaches a scaled one-soliton 
centered at a; = 0. 

By contrast, in the numerical simulations of moving solitons, the defect mode 
is initially forced by the tail at the soliton's leading edge, for which the soliton's 
internal frequency is important. The defect mode may grow by resonantly 
extracting energy from the tail, before the bulk of the soliton even reaches the 
defect. Consequently, when the soliton reaches the defect, the defect mode is 
large enough that the two modes may interact, and energy may flow from one 
mode to the other. In the next section, we introduce an ordinary differential 
equation model of this interaction. 



5 A model of soliton-defect interactions 

Forinash, Peyrard, and Malomed [1] studied the interaction of a soliton with 
a linear defect mode using a collective coordinate ansatz to derive a set of 
approximate equations to describe the evolution of a finite set of variables 
that characterize the two modes. Their model yielded a complicated set of 
differential-algebraic equations which was difficult to understand analytically. 
Here we modify and slightly simplify their ansatz. In particular, we approx- 
imate the solution, u, as the sum of a time-modulated soliton and a time- 
modulated nonlinear defect state: 

u ^ us{x;r],Z,V,4>) + UB{x;a,4>,ip), (5.1) 

where us is a generalized soliton of the form 

us = r;sech(r/a; - Z) e'^^''-^\ (5.2) 

and uy) is a generalized bound state of the form 

md = asech {a\x\ + tanh"^ (,-'^4>+^\ (5.3) 

In (5.2, 5.3), the variables Z, 0, a, and ip are all allowed to depend on t. 
Note that ||'U5|||2 = 2rj and ||-U£)|||2 = 2(a — 7). The Lagrangian functional of 
the NLS equation is evaluated at u (by integrating over the spatial domain), 
and the resulting function is then interpreted as an effective finite-dimensional 
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Lagrangian whose Euler- Lagrange equations determine the evolution of these 
variables. Instead of seeking stationary configuations of the (true) Lagrangian 
with respect to admissible variations of the field variable u{x,t), we consider 
only variations of the type allowed by the time-dependent variables of the 
ansatz. This technique, often known as the collective coordinate or variational 
method, has a long history and is well summarized in a recent survey by 
Malomed [19]. 

The NLS Lagrangian (cf. (3.2)) is 

/oo ^ 11 
-{u*ut - uul) - - \u^\^ + -\u\^ + -i5{x)\u\^dx. 
-oo z z z 

We evaluate this integral for our ansatz (5.1-5.3). Unlike [1] and most other 
analyses of this type of which we are aware, we do not assume that a is small, 
and we include high-order terms involving a. We do make the simplifying 
assumption that all interaction between the modes its ^-nd u-o takes place 
through terms involving the delta function. This is, in part, justified because 
the other interaction terms involve oscillatory integrals which will average 
out to be much smaller than the terms retained. ^ The resulting effective 
Lagrangian is given by 

Aff = 2?70 - 2ZV + 2(a - 7) (0 + V;) + \rf - r)V^ + + ^ sech^ Z 

-\- 2^r]J a? — 7^ sech Z cos ■0. (5.4) 



This Lagrangian has an associated Hamiltonian, conserved by the Euler- 
Lagrange (and Hamilton's) equations: 



H = — rf + riV'^ - '^rf sech^ Z - 2-^7] \J a? - 7^ sech Zcosil). (5.5) 

3 3 

Since the Hamiltonian is independent of 0, its conjugate momentum 

= ^ = 27? + 2(« - 7) (5.6) 

0(f) 

is conserved by Noether's theorem [20]. Notice that =||'^5|li2 + ||'?^D|li2- This 
conservation law is the analogue of the conservation law for NLS, (3.3). 
The phase space of this three degree-of-freedom system can be therefore be 
expressed as the cross product of the reduced 4-dimensional (Z, V, a, ■0)-phase 

^ Here we are essentially anticipating that this effective Lagrangian is a normal form 

to which the 'exact' one is equivalent up to change of variables. The assumption 
that the integrals are oscillatory is violated when \V\ ^ 1, i.e. when the soliton 
is stalled. The assumption that the overlap integral is small will be violated when 
|Z| <C 1, i.e. when the soliton is close to the defect. 
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space and the 2-dimensional (77, 0)-phasc space, with trivial dynamics on 
the latter, determined from the reduced system via (5.6) and a quadrature 
(cf. (5.9-5.10) below). In the analysis that follows, we may therefore regard 
the constant of motion 

c^rj(t) + a(t) (>7) (5.7) 

as a parameter determined by the initial conditions, and study evolution on 

the reduced {Z,V,ip,a)-spacc. Here and henceforth c denotes this constant, 
and should not be confused with its conventional usage to denote the speed 
of light. 

We note further that the resulting equations arc invariant under the rescaling 
a H- > 7a, c H- >■ 7c, t 1-^ t/'j'^: implying that, without loss of generality, the 
parameter 7 may be set to equal to one. This symmetry belongs only to 
the reduced system of ODEs, not to the original NLS equation (2.1), but 
motivates our decision to perform simulations only with 7=1. Then c is 
the only parameter remaining in the evolution equations, and these equations 
are in fact canonically Hamiltonian for the 'scaled' Hamiltonian H/2, with H 
of (5.5) written in the Lagrangian coordinates: 

C3 



H = h (c — a) (ca + V'^ — {c — a) sech^ — 2(c — a)Va^ — 1 sech Z cosip. 

^ (5.8) 
The final reduced equations, to be studied in Section 7 below, are now: 



Z={c-a)V, (5.9a) 

V — — (c — a)^ sech^ Z tanh Z — (c — a) vo^ — 1 sech Z tanh ZcosV', (5.9b) 

— 2ca — V"^ . , , 2 (2a^ — ca — 1) , ^ , . ^ ^ 

-0 = h (c — a) sech Z H — sech Z cos ■0, (5.9c) 

2 Vo^ — 1 

a = — (c — 0)^/0^ — 1 sech Z sin ■0, (5.9d) 
with (p evolving according to: 

(j) — {c — a) sech^ Z — V — 1 sech Z cos -0. (5.10) 



Before analyzing these ODEs, we describe numerical experiments that reveal 
interesting interactions between the soliton and the defect mode, and suggest 
specific questions. 



6 ODE Simulations of the initial value problem 



We now describe a set of numerical experiments for the ODE (5.9) analogous to 
those discussed in Section 4 for the PDE. We initialize a soliton at Z(0) = —Zq 
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with \Z()\ ^ 1, set equal to 20 for these simulations. The velocity parameter 
is set to V^(0) = Vin, propagating rightward toward the defect. For subsequent 
comparison with the PDE simulations, we note that the variable V in (5.9) 
is related to the soliton velocity v of (2.2) via v — rjV by (5.9a) and (5.7). 
The soliton amplitude is set to 7^(0) = rjo, and we assume there is no energy 
initially in the defect mode, so a(0) =7 = 1. Thus the relation (5.7) fixes the 
constant c. (Because of the singularity of Equation (5.9c), we set a(0) = 1+e, 
where £ -C 1; values between e = 10~^ and s = 10~^ were used in these 
computations. For some of the effects seen, e needs to be set to the small end 
of this range to get stable results, due to the singularity in (5.9c). ) Finally, we 
set ■0(0) = 0. Note that, for large \Z\, the ip dependence becomes exponentially 
weak, so we do not expect strong '^-dependent effects. 

We choose a representative set of ?7(0), or equivalently, c values. For each fixed 
7^(0) we allow T^n to vary over a range of values, and numerically integrate (5.9) 
until the soliton center, Z{t), reaches the defect {Z{t) ~ 0), (eventually) exits 
the defect region, and reaches Z = +Zo or Z = —Zq at, say t = T. By a variant 
of the Poincare recurrence theorem, as in [5], we can show that the soliton 
must eventually escape any bounded set containing the defect; cf. Section 7.1 
below. We then plot the soliton's outgoing velocity parameter V^ut (related to 
the physical velocity by (5.9a)), and the amplitude (a(T) — 1) of the defect 
mode as functions of Vin. For Zq sufficiently large and e sufficiently small, the 
initial value of the phase difference ip was indeed found to be unimportant in 
determining the values of V{T) and (a(T) — 1). In interpreting these results, 
it is important to realize that the dynamics takes place in four-dimensional 
phase space, and that the figures merely show projections of trajectories on 
lower dimensional subspaces. 



6.1 Experiment 1: large r]. 

We observe several distinct types of behavior. The behavior for 77 = 4 is 
shown in Figure 6.1. In this case we find a sharp change in behavior at a 
critical velocity parameter Vc ~ 0.55. Above this velocity, sohtons pass through 
the defect without significant interaction, merely decreasing their velocities 
and transferring a little energy to the defect mode. Below Vc, however, the 
soliton interacts with the defect, oscillating within the defect region a finite 
and apparently random number of times before being ejected either to the right 
or the left. It is also striking that the amount of energy remaining in the defect 
mode seems to be restricted to two levels: either a — iRiSora— l<Sl. This 
behavior is apparently governed by the structure of the invariant manifolds of 
degenerate fixed points at \Z\ = 00. Since there exist solutions to (5.9) with 
Z bounded, as well as solutions which approach Z = 00 with V > 0, these 
must be separated by solutions for which Z — > 00 while y — > 0. Along with 
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Fig. 6.1. (Top) Outgoing velocity vs. incoming velocity of solitons incident on defect 
via ODE simulation with rj = A. (Bottom) Amplitude a — 1 of nonlinear defect mode 
after passing of the soliton. 



the apparantly arbitrarily fine structure of transmission and refiection zones, 
Figure 6.1 further shows wide reflection windows of the type reported in many 
previous studies, e.g. [21-24]. 

In Section 7.1 we will investigate the stable and unstable manifolds that are 
responsible for this behavior, but we give a brief preview here. Figure 6.2 shows 
the {Z, V) projections of trajectories with nearby initial velocities on either 
side of Vc {V{0) — 14=1= -002). The solid curve comes close to an orbit apparently 
asymptotic to {Z, V) = (oo, 0), and then turns back and is transiently captured 
before eventual reflection, while the dashed curve approaches infinity with 
V bounded above zero, and is transmitted without further interaction. This 
'separatrix' behavior is repeated in the multiple transmission and refiection 
windows for V < Vc (Figure 6.1) and is reminiscent of that found in our 
earlier study of an ODE model of kink-defect interaction in the sine-Gordon 
equation [5], and shown there to be related to the homoclinic tangle formed 
by transverse intersection of stable and unstable manifolds of periodic orbits 
at \Z\ = oo: compare Figure 6.1 with Figure 3.2 of [5]. 

The ODEs' behavior should be compared with the direct numerical simula- 
tions of solitons with = 4, as shown in Figure 4.1. The critical velocity 

is overestimated by 24% {vc = VjVc = 2.21 cf. 1.78 for the PDE), and the 
PDE simulations show no evidence of the fine structure of transmission and 
refiection zones below Vc- This difference, partially due to neglect of radiation 
damping in the ODEs, is discussed in section 8. 
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Fig. 6.2. (Z, V) phase space plots of captured and transmitted trajectories just above 
(dashed) and below (solid) critical velocity. 
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Fig. 6.3. (Top) Outgoing velocity parameter V vs. incoming V of solitons incident 
on defect via ODE simulation with 77 = 2. (Bottom) Amplitude a — 1 of nonlinear 
defect mode after passing of the soliton. 

6.2 Experiment 2: medium 77. 

In Figure 6.3, computed for 77 = 2, we see a quite different picture. Here and 
in the third numerical experiment, below, —rf/'2. < —7^/2 = —1/2 and reso- 
nant interactions can and do take place (cf. Figure 2.1). In this case there is 
no transition between transmission and reflection: the sohton travels mono- 
tonically rightward and is always transmitted without transient capture or 
oscillations about the defect. More strikingly, the output velocity appears to 
approach a finite limit V^ut ~ 1-17 ^ 0, while the amount of energy 

captured approaches a — 1 ~ 1.26 (all the soliton's energy would be captured if 
a — 1 = 2). As the initial velocity increases, so does the output velocity, while 
the energy transferred from the soliton to the defect mode decreases. This is 
not surprising, since the duration of the interaction decreases with increasing 
soliton speed. 
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Fig. 6.4. (Top) Outgoing velocity vs. incoming velocity of solitons incident on the 
defect via ODE simulation with rj = ^. (Bottom) Amplitude a — 1 of nonlinear 
defect mode after passing of the soliton; note vertical scale is O(10~^). 



6.3 Experiment 3: small rj. 



For Figure 6.4, we set = 0.5 < 7; in this case the sohton is reflected if T^n 
hes below a critical velocity Vc ~ 0.51 and transmitted if V^n > We also see 
that the defect mode has a final amplitude of order 10^^, absorbing little of 
the soliton's energy. We note that —7^/2 = — 1/2 < —rf/2 corresponds to the 
region in Figure 2.1 in which the soliton has no resonant defect mode 'partner' 
with the same temporal frequency, and hence that appreciable interactions are 
unhkely [2]. 

Figure 6.5 shows evidence that the solution passes near Z = V ~ presum- 
ably approaching and leaving the neighborhood of a hyperbolic invariant set 
on this subspace. In Section 7 we shall show that this subspace indeed contains 
a fixed point of saddle-center type surrounded by a family of periodic orbits 
whose stable manifolds serve as separatrices. In Figure 6.5 we also show pro- 
jections onto (Z, l/)-space of the numerically determined stable and unstable 
manifolds of the fixed point, along with two trajectories: one with asymp- 
totic velocity larger than the limiting value on the unstable manifold, which 
is transmitted, and one with asymptotic velocity smaller than the limiting 
velocity, which is reflected. 
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Fig. 6.5. Phase space projections of trajectories above (solid) and below Vc (dashed), 
showing strong evidence for a saddle point. Projections of stable and unstable man- 
ifolds of the saddle on Vq arc shown as dotted curves. 

7 Analysis of ordinary differential equations 

The numerical experiments described above reveal two broad types of behav- 
ior. For small values of 77, the soliton traveling near the critical velocity appears 
to approach a hyperbolic fixed point or periodic orbit on {Z, V) = (0, 0): see 
Figure 6.5. For larger values of 77, it oscillates about {Z, V) = (0, 0) as if around 
an elhptic fixed point, and also follows orbits asymptotic to (Z, V) — (±00, 0) 
before turning: sec Figure 6.2. To interpret these observations we now analyze 
the ODE system (5.9), seeking a (partial) understanding of its global phase 
space structure. 

We first note that the sets a = c and a = 1 are invariant for the flow (although 
the vector fleld is singular on the latter), and bound the physically admissible 
region. When a = 1 all the energy resides in the soliton; when a = c it all 
resides in the defect mode. Letting I = [l,c] denote the closed interval, the 
phase space of Equation (5.9) is (Z, V, t/j, a) e x 5"^ x X. We also note the 
following (reversibility) symmetry group under which (5.9) is equivariant: 

{Z, V, ip, a, t) {-Z, V, 271 - i/j, a, -t), (7.1a) 
{Z, V, ^, a, t) {Z, -V, 2% - ^p, a, -t). (7.1b) 

We shall use this below. 

There is a family of solutions at Z = ±00 with V, a, and ip constant, which 
correspond to the uncoupled propagation and oscillation of the two modes 
when the soliton is infinitely far from the defect. The subset of these solu- 
tions with V = form a degenerate family of periodic orbits 'at infinity' 
parameterized by a = Ooo and filling the annulus (or finite cylinder) 

Poo = {(V',a)l^ = 0, |Z| = oo}. (7.2) 

We note that by (5.9c) and (5.7), = (r^^— a^)/2 on Voo, so that the frequency 
of these orbits is nonzero provided a^o 7^ V or, equivalently, aoo 7^ c/2. As in [5] 
we may employ a transformation of the form q — sech Z, p — V to bring these 
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orbits to the origin in {p, g)-space, and then apply McGehee's stable manifold 
theorem [25] to prove the existence of invariant manifolds for the fixed point 
(P) q) = (0, 0) in an appropriate (local) Poincare map. This shows that each 
periodic orbit in Vqo has two-dimensional stable and unstable manifolds, so 
that W^{Voo) - the stable manifold of Voo itself - is three dimensional and 
hence locally separates the four-dimcnsional phase space. Indeed, W^lVoo) 
separates orbits that escape to infinity (transmitted solitons) from those that 
are refiected to interact with the defect mode again. 

Studying the analogous sine-Gordon kink-trapping problem in [5], we used 
isoenergetic reduction and Melnikov's method [26,27] to prove that the sta- 
ble and unstable manifolds of each periodic orbit, restricted to their com- 
mon energy manifold, intersect transversely, and hence that Smale horseshoes 
exist [27]. We then appealed to phase space transport theory [28,29] to un- 
ravel the structure of sets of initial data that are transiently captured before 
eventually being transmitted or reflected. We proceed in the same manner 
in Section 7.1, although we have to introduce an artificial small parameter, 
and the set of solutions to which the standard reduction procedure applies is 
limited, since it requires that the frequency ip not change sign (in the process 
one replaces time by and this holds only for large c, depending on a; cf. 
Equation (5.9c). In particular, it does not hold for many physically relevant 
parameter values, including those corresponding to initial data with (almost) 
all the energy in the soliton (a ~ 1). Nonetheless, the analysis does provide 
some understanding of the large rj simulations of Section 6.1. 

There are further invariant manifolds that play an important role in the fate 
of solutions. They belong to orbits on a second annulus 

Vo^{i^,a)\V^O, Z = 0} (7.3) 

that is also invariant under the fiow, and on which the ODEs are integrable. 
Solutions on Vq correspond to a soliton stalled over the defect and periodically 
exchanging energy with the defect mode. The orbit structure on Vq depends 
upon c and may be derived from the level sets of the Hamiltonian function H 
restricted to Vq: 

H-p^{ip, a) = -— + {c-a) [ca - (c - a) - 2 Va^ - 1 cos ip) . (7.4) 

As noted above, the boundaries a — c and a = 1 of T'o are invariant, and the 
flow is singular on the latter, which contains two degenerate saddle points at 
t/j = 7r/2, 3n/2. There is a unique fixed point (vr, a*) on ip = n surrounded by 
periodic orbits which limit on heteroclinic orbit (s). As shown in the first section 
of the Appendix, for all c (> 1), (tt, a*) is a saddle-center, with positive and 
negative real eigenvalues whose eigenvectors point out of Vo- For c < 2.214 .. . 
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Fig. 7.1. Level sets of Hq in the invariant plane "Pq for (a): c = 5 (r/ = 4); (b) 
c = 3 (rj = 2) and (c) c = 1.5 (rj = 0.5). The dashed curve indicates the boundary 
of the orbits accessible from infinity, as described in the text. 



this is the only equilibrium; for c > 2.214 . . . two further fixed points, a center- 
center and a saddle-center, appear on ^0 = 0; restricted to Vo these are a 
center and a saddle, whose separatrices interact with the stable and unstable 
manifolds of the degenerate saddles on a = 1 in a heteroclinic bifurcation [27] 
at c ^ 3.21 as c continues to increase. Ho takes its minimum value — c"^/3 on 
c = a, its maximum at (tt, a*), and the value —c^/3 + c— lona = l and the 
invariant manifolds emanating from it. Figure 7.1 shows these distinct cases. 
On this figure we also show the level set with Hamiltonian value equal to that 
of a 'pure' soliton stalled at infinity: H{\Z\ = oo, V = 0,a = = — c^/3 + 
c(c— 1). Since any incoming soliton with nonzero speed has H > — c'^/3+c(c— 1) 
(see (5.8)), this curve bounds the set of accessible orbits on Vq: a disk centered 
on (tt, a*). 
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7.1 Stable and unstable manifolds of Voo 

We first observe that we may define a local three-dimensional cross section [27] 
for the fiowof (5.9): 

= {{Z, V, i^, a)\a e (a*, c),i^ ^ tt}. (7.5) 

We verify that the flow is transverse to E^r in the second section of the Ap- 
pendix. 

Since the Hamiltonian (5.8) is time-independent, its value is conserved by 
solutions of (5.9), which are therefore constrained to lie on three-dimensional 
surfaces 

H(y, I]c) — h — const, 

determined by the initial data. The variables p and /, defined below, denote 

coordinates which are conjugate to coordinates V and ip. As shown in the 
final section of the Appendix, this permits a further reduction in dimension. 
Specifically, since the cross section intersects level sets of H transversely, 
E^ may be used as a two-dimensional cross section for the fiow restricted to 
constant H surfaces. It is on this cross section that we will portray the stable 
and unstable manifolds W^''^{Voo)- 

To approximate these manifolds we first add an artificial coupling parameter 
(1 to the Hamiltonian of Equation (5.8): 

H = h (c — a) {ca -\-V'^ — {c — a) sech^ — /x2(c — a)\/a^ — 1 sech Z cos if) . 

^ (7.6) 
For the case at hand, /i. = 1, but we shall assume /x -C 1 and perform a 
perturbative analysis, subsequently appealing to continuation to extend to 
^ — 1. The variables V and ■0 are canonical 'positions' for this Hamiltonian 
with conjugate 'momenta' p = —2Z and / = 2(a — 7) = 2(a — 1), and / 
and ip are action- angle variables for the 'second' degree of freedom. In these 
canonical variables, the Hamiltonian (7.6) assumes the form: 

H^H,{V,p■,I)+,,H^{p■,I,^p). (7.7) 

The formal discussion of the Melnikov integral will refer to these canonical 
variables, although in both the computations to follow, and in numerical sim- 
ulations, it is more convenient to work with the original variables {Z, V, ip, a). 

For II — 0, the 'unperturbed' Hamiltonian Hq is independent of ip, and the 
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ODEs reduce to 



Z={c-a)V, (7.8a) 

V ^-{c-a)'^sech^ZtanhZ (7.8b) 
— 2ca — V"^ 

h(c-a)sech^Z (7.8c) 

d = 0. (7.8d) 



The position and velocity Z and V evolve as a particle in a potential well, 
with the strength of the well dependent on a — Oq, which is unchanging. 
The solution set comprises bounded periodic orbits, and unbounded orbits 
where \Z\ ^ oo with finite speed. In between is a separatrix, on which 
Z = sinh"^ Vo^'^t^ where r/o = c — Oq is the amplitude of the soliton. In the 
unperturbed dynamics, the angular displacement, determined by integration 
of (7.8c) along the separatrix, is given by 

^ = ^0 + + ^ arctan rj',/'. (7.9) 

The separatrices are homochnic orbits to a periodic orbit P{t) e Voo with 
a = Oq and {Z, V) = (±oo, 0). Note that unlike in [5] and related problems in 
celestial mechanics (e.g. [30]), in this case the generalized 'position' (soliton 
speed, V) goes to zero and the conjugate 'momentum' (soliton position, Z) 
goes to infinity. 

We ask if any members of this continuum of homoclinic orbits persist when 
H ^ 0. Now if W^{f3) crosses the p = axis at the point {p = 0,V,ip = 0,a,t = 
0), then, by the symmetries (7.1), W'^{P) must also pass through this point. 
To establish existence of a (perturbed) homochnic orbit it therefore suffices to 
show that either W^{P) or W"^{l3) intersects p = 0. However, to demonstrate 
transverse intersections requires a more delicate calculation, appealing to per- 
turbation theory for small n and a result due to Melnikov, the derivation of 
which is outlined in the final section of the Appendix. Specifically, we have: 

Theorem 2 Let h > and Vlo{t) = Let {Ho,^}denote the Pois- 

son bracket^ of Ho{V^,p°) and iJi(V,/, /°)/Oo(V^°,P°) evaluated along 
V^{t) and p^ it). Define the Melnikov function 

M(^o) = /_^{i/o,^}(V^,P°,^°,/°)(t) dt, (7.10) 

and assume that M{iPq) has a simple zero and ilo(^) 7^ 0. Then for /i > 
sufficiently small, the Hamiltonian system has transverse homoclinic orbits on 
the energy surface H — hP. 

2 r F r^i def OF dG dF dG 

i^^'^s — dv dp dp av 
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We note that the usual reduction process and Melnikov integral are meaningful 
only as long as i/j is monotonic with respect to t {ip = Q ^ 0), so that the 
global cross section S^g referred to in the Appendix may be defined. However, 
Holmes [31] has extended this analysis to the case where ip is allowed to change 
sign in a bounded region in the 'middle' of the unperturbed homoclinic orbit, 
but under the condition that ip be monotonic sufficiently close to the periodic 
orbit (3. Direct substitution of the unperturbed orbit Z — sinh r)Q t into 
Equation (7.8c) yields 

By our ansatz rjo > 0, and so the condition that Qq not change sign throughout 
is 

(1 <) ao < 770 or ao > ^/vo+Vo; 

however, if we appeal to [31], we need only exclude a neighborhood of the 
degenerate set 

ao = 770 

to ensure that fio 7^ near /3. 

We now sketch the computation needed to verify the remaining hypothesis of 
Theorem 2. Since Hi is F-independent, the Poisson bracket reduces to 

^ Hi] 1 dHodHi Hi fOHoOno dHodVlo 

Ho, 



' ^0 i (v,p) ^0 dV dp fll \ dV dp dp dV J ' 

After computing partial derivatives (cf. (7.8)) and substitution of the unper- 
turbed separatrix solution in the Hamiltonian of (7.6), some manipulations, 
and appeal to odd- and evenness properties, the Melnikov integral may be 
written as: 

MiiPo)=4M^lir (-(^o-"o)^(^)+^°^"^(^))f^^'^^tlsm^o, 

(7.12) 

where 

n2 ^2 



g{t) = Jl + rilf^ and Q{t) = — — —t + arctanr/l^- 

2 



This integral cannot, in general, be computed explicitly, unless ^/fjo is ratio- 
nal. However, since the integrand is an analytic function of t with no essential 
singularities, it can have only isolated zeros. Therefore, except for special val- 
ues of (r7o,ao), Ai{iljQ) has only simple zeros where ipo = mr, and Theorem 2 
implies that, for small values of /x, there exist transverse homoclinic orbits to 
infinity. 



26 



Fig. 7.2. Stable and unstable manifolds of periodic orbits /? G Too, shown via the 
Poincare map defined on the cross section S^r. (a) fx = 0.25; (b) /x = 0.5, corre- 
sponding to the energy level for = 4 and V = 0.4. 

Figure 7.2 shows numerical computations of the stable and unstable manifolds 
of periodic orbits in Too for two values of ii in the same energy surface H — 
—20.69, corresponding to the energy of the system with a soliton with rj = 
4 and V = 0.4 starting at \Z\ = oo, with a = 0. They are illustrated as 
curves lying in the cross section E^r introduced at the beginning of this section. 
(Since each periodic orbit /3 in Voo is a one-dimensional circle, W^{P) and 
are each two-dimensional, and so intersect suitable cross sections to 
the flow in one- dimensional curves.) The transverse intersections are clear from 
the figure. As described in [5], phase space transport theory [28,29] may be 
used to analyse the capture, transmission and reflection dynamics implied by 
transveral homoclinic points such as those of Figure 7.2. Here we provide a 
brief review; for a more complete explanation, see [5]. 

We consider W^{Z = —oo), the unstable manifold of Z = —oo, and W'^{Z ~ 
oo), the stable manifold of Z = oo, which intersect transversely in a point 
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in the top center of the figure. (Both cases show such an intersection, although 
the phenomenon is clearer in the lower figure.) The union of the point the 
portion of W^{Z = — oo) to the left of Qq, and the portion of W^{Z = oo) 
to the right of go form a boundary between the upper and middle regions of 
the plane. A similar boundary exists in the lower half plane. Stretching off to 
the left from the point go is a series of lobes lying in the upper region, each of 
which is the image under the map P^^ of the lobe to its left. Between each such 
pair of lobes in the upper region, and below W^{Z = — oo), lies a lobe in the 
middle region. Counting both sets of lobes, the image of any given lobe under 
is the second lobe to its right. In particular, the image of the nearest upper 
region lobe to the left of go is a lobe located in the middle region. Similarly, 
the image of the middle region lobe immediately to the left of go lies in the 
upper region. This nearest upper lobe and the neighboring middle lobe form a 
turnstile through which phase space points are transported from the upper to 
the middle region, and from the middle to the upper region. A similar turnstile 
in the lower half plane transports phase space between the lower and middle 
regions. 



Trapping takes place when an initial condition that lies in the sequence of lobes 
in the upper-left quadrant is mapped from the upper region to the middle 
region. Refiection or transmission occurs because eventually, as the interior 
lobes are successively stretched and folded, the image of this point will, with 
probability one, lie in a turnstile exit lobe (area preservation of the symplectic 
map Ptt guarantees that no open set contained in a preimage of an incoming 
turnstile lobe can be trapped for all future time [5, Proposition 1]). If the 
point's image exits into the upper region, it is transmitted. If it goes into 
the lower region, it is refiected. This may be seen as an analogy with phase 
space transport in a two-bend horseshoe map, as is shown in [5] . The iterated 
preimages of the turnstile lobes forms a fractal structure that is responsible 
for the multiple refiection and transmission windows of Figure 6.1. 



7.2 Stability of orbits on Vq 



To determine the stable and unstable manifolds of Vq we must first determine 

the stability types of orbits within it, to perturbations out of Vq- By continuity, 
the periodic orbits immediately surrounding the saddle center (vr, a*) are also 
of saddle type with respect to such perturbations, but the stability types of 
other periodic orbits must be determined via Floquet theory [32]. 
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On Vo the ODEs reduce to: 




— ca + {c — a) + 



(2a2 -ac-1) 



COS'0, 



(7.13a) 



a — — 



(c — a)y/a'^ — 1 smijj. 



(7.13b) 



Typical phase portraits of (7.13) are shown in Figure 7.1 above. 

Consider the solution S* to the initial value problem of a soliton starting from 
|Z| = oo with finite velocity and zero energy in the defect mode, i.e. a — 1 
and c = 77 + 1 ^ . <S* is confined to the level set 



of the conserved Hamiltonian and hence, if it approaches Vo, can only interact 
with orbits having the same iiT-value. In particular, since the maximum H 
value for orbits on Vq is assumed by the fixed point (tt, a*), there is a critical 
velocity V^*^^ above which the solutions S* have more 'energy' than any orbits 
contained in Vq, and thus must remain bounded away from it. Similarly, the 
minimal value H = — + c(c — 1) of orbits S*, assumed when = 0, 
bounds the set of accessible orbits on T^o, as shown by the dashed curve on 
Figure 7.1. 

Consequently, for each V^o e [0, V^^), we find a periodic orbit Sq G Vq with 
the same Hamiltonian value H = as S* and determine its stability by exam- 
ining the hnearization of the full system (5.9) about Sq = {0,0, tpp(t),a-p{t)). 
The stability of such an orbit is given by the eigenvalues of the monodromy 
matrix: the fundamental solution matrix of the linearized differential equa- 
tion, evaluated at one period of oscillation. Let <So = {Z, V, ip, d) solve this 
linearized ODE, which is block-diagonal, with the {Z, V) components decou- 
pling from the (ip, d) components. The cigenspace of the latter coincides with 
Vo and hence they have eigenvalues of unit modulus, as one expects from the 
integrable structure of Figure 7.1 (in fact A = 1 with multiplicity two and a 
single eigenvector). 

Perturbations perpendicular to Vq satisfy 



^ Practically, as in Section 6, the orbit is started at some \Zo\ 3> 1 with a = 
1 + e, £ < 1. 



^ = -y + (c-l)(c + = ^c 



00 



(7.14) 





{c-ar{t))V, 



(7.15) 
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Since this may be written in the form 




fZ 
A2,i{t) j V 



) 



the Floquet theory for Hill's equation is applicable [32]. The product of the 
Floquet eigenvalues must be one, and their sum is given by the Floquet dis- 
criminant. If this is greater than two in absolute value, the periodic orbit is 
hyperbolic, i.e. unstable; if less than two, the orbit is elliptic, i.e. neutrally 
stable, transverse to Vq. We approximate these discriminants by first numeri- 
cally integrating the orbit {a'p{t),ip-p{t)) for one period, and then computing 
the fundamental solution matrix for the linearized system (7.15), using inter- 
polated data for the coefficients. 

Since each periodic orbit in Vq corresponds, via its Hamiltonian level, to a 
velocity at |Z| = oo, we plot in Figure 7.3 the Floquet discriminants as 
functions of Voo for the three examples of Section 6.3. In the first case {rj — 4) 
there are two regions each of stability and instability, and the velocity range 
shown in Figure 6.1 corresponds to periodic orbits of elliptic type, consistent 
with the intuition from Figure 6.2, that solutions oscillate about Vq. The 
critical velocity dividing capture and transmission, identified in Figure 6.1 
is indicated by an asterisk. The second case oi r] — 2 encompasses a range 
of stability and one of instability, and in the third (and simplest) case, with 
1] = ^, all orbits are hyperbolic. It is notable that V^^^^^^ = 0.51 in this case, 
approximately equal to the critical velocity for transmission and reflection 
seen in Figure 6.4 (but see the detailed analysis in Section 7.3, below). 

Some interesting implications can immediately be drawn from the stability 
types of orbits in Vq implicit in Figure 7.3. Recall that these periodic orbits 
correspond to a state in which a soliton, stalled over the defect ai Z — 0, 
exchanges energy periodically with the defect mode. Hence, if soliton and 
defect parameters are chosen consistent with a stable region (eg., below T4o ~ 
0.9 in case rj = 2 and for Voo < 1 and Voo G (1-4, 1.8) in case 7] = 4), and 
the soliton is initialized at the defect or somehow introduced into it, perhaps 
by temporarily destabilizing the relevant orbit, it will remain trapped under 
small perturbations. Such stable trapped states do not exist for smaller rj. 

In the case rj = 0.5, the initial condition S* = (— oo, Vc = 0.501, l,ipo) hes on 
the same Hamiltonian surface as the fixed point (7r,a*). The stable and un- 
stable manifolds of this fixed point are only one dimensional and thus cannot 
separate refiected from transmitted orbits. However, the stable and unstable 
manifolds of (the accessible disk on) Vq are three dimensional, and are conse- 
quently able to divide the phase space into disjoint regions. We must therefore 
compute the stable and unstable manifolds of the accessible periodic orbits on 
Vo as well as of the saddle-center (tt, a*). 
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Fig. 7.3. The Floquet discriminant vs. incoming velocity for 77 = 4, 2, 2 respectively 
(top to bottom). In the top panel, the critical velocity separating capture from 
transmission is marked with an asterisk. 

1.3 Stable and unstable manifolds ofVo 

Appealing to the symmetries of (7.1) and the fact that the orbits of interest 
are reflection-symmetric about ip — tt, we need only compute one of the four 
branches of W'^'^{So) for each of the saddle type periodic orbits Sq G Vq- To do 
this we first compute each periodic orbit So starting at a point (a, tp) = (oq, tt) 
where gq > a*, the saddle-center. We interpolate these with 64 equally-spaced 
points (with respect to time). At each of these 64 points (ao,'0o), the funda- 
mental solution matrix is computed as in Subsection 7.2. Fourier interpolation 
is used to compute the coefficients {a-p,ip-p) at intermediate times, so the or- 
bit So need only be computed once. At each point on the periodic orbit, we 
compute the unstable eigenvector of the monodromy matrix, vo = {Zq, Vq)^. 
We normalize it so that \vo\ = 10~^ and solve the full system (5.9) of ODE's 
with initial conditions (Zq, Vq, Oq, V'o); stopping when \Z\ — 20. 

Let W^{Vo) = UhoW^{So) denote the set of stable manifolds of the accessi- 
ble hyperbolic orbits in Vq. W''{Vo) is three dimensional, and is locally (near 
\Z\ = 0) fohated by two-dimensional cylinders, each of which is a local sta- 
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ble manifold of some <So. We therefore expect W^{Vo) to intersect the three- 
dimensional cross sections of initial data Sj. = {{V. iL\a)\Z = ±Zq, \Zo\ ^ 1} 
in two-dimensional sets, which should in turn separate sets of initial data 
giving rise to solutions that pass the defect from those reflected by it. 

Figure 7.4 shows the results of computations for the third of the three cases of 
Section 6: t] = 0.5, and for a further case, with slightly larger r] = 0.75. Note 
that the sets Ej. ^ M. x x 2 are periodic in ijj. As might be expected from 
the experiment of Section 6.3, near a — 1 the surface W^{So) n E_ is a graph 
over the a) annulus: all orbits starting at points below it are reflected, 
and points starting above it are transmitted. Further from a = 1 the surface 
develops folds; these become more pronounced for higher rj (c), as in the lower 
panel of Figure 7.4. The surface describes the critical velocity as a function 
of phase t/j and amplitude a — 1 of the defect mode. Note the weak phase 
dependence, particularly as a — > 1, and that the surface approaches a = 1 at 
V ~ 0.51 in the case r] = 0.5 (upper panel); this is the critical velocity found 
in Section 6.3. Initial data on this surface corresponds to trapping (recall that 
the accessible orbits in Vq correspond to solitons pinned at the defect). 

In interpreting this figure it is helpful to note the following facts. Individual 
two-dimensional components W^{Sq) ofW^iVo) intersect sections at Z = ±Zo 
for small \Zo\ in homotopically trivial (contractible) circles, but as \Z\ (and 
the time of flight) increases, particular solutions belonging to W^{So) can 
pass arbitrarily close to the degenerate saddles on a = 1 at = 7r/2,37r/2 
(cf. Figure 7.1 (bottom)). This effectively separates neighboring solutions and 
stretches their phase (ip) angles over a range exceeding 27i. The result is that 
the corresponding sets W''^{So) H S_ extend around the component in a 
homotopically nontrivial manner. Only those components ofW'^{Vo) very close 
to W^{7r,a*) remain contractible; these can be seen in Figure 7.4 near tjj ~ 
7r,a — 1. 

The 'outer' (lower V, higher a) boundary of the computed portion of W''^{So) H 
E_ is limited by numerical issues: it is impossible to compute with uniform 
accuracy as velocities approach zero, since the time of flight grows without 
bound; however, it appears that velocities do decrease to zero as a increases. 
For example, numerical experiments like those of Section 6.3 indicate that all 
orbits launched with positive velocities, no matter how small, and a > 1.1, are 
transmitted. The surface therefore intersects V = 0. 

We have been unable to reliably compute invariant manifolds of Vq for the 
medium and large r) cases. Preliminary studies suggest that, as c increases, the 

set W''^{Vo) n E_ 'separates' from the plane a = 1, so that nearby initial data 
all lie in the transmission zone (cf. Figure 6.3, which indicates that all solitons 
are transmitted for rj = 2, regardless of their initial velocities). However, as 
c continues to increase, the stable manifold W^{Voo) evidently invades E_, 
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leading to the complex behavior of Figure 6.1. In particular, the increased 
folding of W^{Vo) n S_ as (or c) increases evident in Figure 7.4 is consistent 
with the existence of a fine (fractal) structure suggested by Figure 6.1. Since 

W''^{Voo) cannot intersect W'^{Vo), we conjecture that, as W'^{Voo) invades S_, 
it must 'align' with the latter (folded) manifold, producing (infinitely) many 
regions of transmission and reflection on any vertical line in E_ above the 
(a, '^)-plane. 
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8 Interpretation and Sumniciry 



In this paper, we have derived a finite dimensional model for soliton-defect 
mode interactions in a nonlinear Schrodinger equation with a point defect. 
Following [1], and allowing for a fully nonlinear defect mode, which by itself 
is an exact solution, we derive a three-degree-of-freedom Hamiltonian system 
that describes the evolution of amplitudes and phases of the soliton and defect 
mode, and the position and velocity of the former. Allowing a nonlinear defect 
mode is important, since it permits resonant energy transfer to occur over a 
range of soliton amplitudes. However, only these two modes are represented; 
in particular, radiation to the continuum is ignored, and multiple solitons are 
disallowed. 

The resulting ODEs may be further reduced to two degrees of freedom, since 
in addition to the Hamiltonian, a second quantity, corresponding to the total 
energy in the two modes, is also conserved. While this system is rather com- 
plex, and indeed is nonintegrable in certain parameter ranges (cf. Section 7.1), 
it possesses two-dimensional invariant subspaces filled with periodic orbits, 
whose stable and unstable manifolds partially determine the global structure 
of solutions. We use this system to investigate the reflection, transmission, 
and trapping of solitons launched 'from infinity,' by the defect, concentrating 
on the case in which the energy initially all resides in the soliton. Numerical 
simulations of the model ODE's reveal three basic types of behavior: 

(1) For large initial soliton intensities, there is a critical velocity above which 
all solitons are transmitted; below this, a complex structure of reflection 
and transmission bands exists, separated by trapping regions that ap- 
parently are of measure zero. The capture and eventual transmission or 
reflection is explained by phase-space transport via turnstile lobes formed 
from parts of the stable and unstable manifolds of orbits 'at inflnity,' cor- 
responding to uncoupled oscillations of the defect mode and a distant 
soliton. 

(2) For moderate initial soliton intensities, all solitons are split into a trans- 
mitted part and a captured defect mode. The transmitted part travels 
to the right monotonically, giving up a fraction of its energy to the de- 
fect mode. The amount of energy transferred to the defect mode is a 
decreasing function of incident velocity. 

(3) For small initial soliton intensities, reflection or transmission occurs for 
almost all initial velocities. Specifically, a unique critical velocity exists 
for each initial phase and defect amplitude below a certain level; this rep- 
resents initial data on the stable manifold of a subset of periodic orbits, 
each of which corresponds to the soliton stalled over the defect, periodi- 
cally exchanging energy with the defect mode. This can be explained by 
the stable and unstable manifold of the manifold Vq, which divide the 
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phase space into two regions. 

We now compare this to the behavior of the numerically computed solutions 
to the original PDE: 

(1) For large initial soliton intensities, there also exists a critical velocity 
separating solitons which are captured from those which pass by the 
defect with little interaction. The ODE prediction of Vc = f]Vc = 2.21 
(with 77 = 4) is in error by some 24% compared to the PDE simulations, 
Vc — 1.78, and unhke solutions to the ODE, captured PDE solitons do not 
eventually escape. The reduced ODE system is Hamiltonian, and thus, 
by a variant of the Poincare recurrence theorem, any solution for which 
Z is unbounded as t — > —00 must also have unbounded Z at a later 
time, with probability 1. In the full PDE, however, radiation damping 
plays a role; there are radiative modes which carry energy away from the 
defect mode. These radiation modes are not incorporated in our Ansatz; 
as in other problems (cf. [5,33]), their inclusion is expected to yield a 
collective coordinate ODE reduction with corrections which can play the 
role of damping. This damping corresponds to energy transfer from the 
soliton-defect mode subsystem to the radiative 'heat bath.' While the 
finite dimensional Hamiltonian ODE reduction leads to trapping for a 
set of data of measure zero, for the reduction perturbed by damping, a 
set of data of positive measure is trapped. 

The ODE model displays two behaviors for solitons below the criti- 
cal velocity. Figure 6.1 reveals the existence of reflection windows: initial 
velocities for which the soliton returns to minus infinity with the same 
intensity and the opposite velocity it started with. Between these reso- 
nance windows arc chaotic regions, where the soliton may oscillate near 
the defect any number of times before being ejected, with apparently ran- 
dom velocity. Reflection windows are a common feature of ODE models 
of this type [23,21] as well as in some PDEs describing soliton-defect 
and soliton-soliton interactions [22,24]. However, the fractal structure of 
the reflection/transmission windows is seen only in the ODE models and 
not in the PDEs. Radiation damping becomes important when the soli- 
ton stays in the neighborhood, and eliminates the chaotic behavior. It is 
shown in [5] that the chaotic behavior can be eliminated by the inclu- 
sion of appropriate radiation terms in the ODE ansatz, which leads to 
damping in the ODE. 

In numerical simulations of NLS soliton-defect interactions, no reflec- 
tion windows have ever been found via PDE simulations. It appears from 
our simulations, and from a form of post-processing of the simulations 
to be described momentarily, that energy may be transferred from the 
soliton to the defect mode, but not vice-versa, so that once a defect mode 
is created, it never gives up its energy to the soliton. In the numerical 
post-processing, we calculate the six ODE parameters of (5.4) by mini- 
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mizing the distance between the numerical PDE solution and the ODE 
ansatz defined by (5.2) and (5.3). In this analysis, for captured solutions 
T] decreases to zero as a grows, so that the soliton is destroyed as the de- 
fect mode is created. This is in contrast with the sine-Gordon simulations 
of [23] in which reflection windows are seen. The difference hes in the fact 
that the two interacting modes in the sine-Gordon experiments are 'topo- 
logically' distinct. In that case the soliton is a kink, which approaches two 
different limits as x ^ ±oo, while the defect mode is exponentially lo- 
cahzed. The kink is defined solely by position and velocity, and has no 
amplitude parameter equivalent to r). When the kink transfers energy to 
the defect mode, it is not destroyed, since it still must satisfy the bound- 
ary conditions at ±oo. This topological constraint forces it to persist, so 
that energy stored in the defect mode may be (re-) converted to kinetic 
energy that pushes the kink away. In the NLS, in contrast, there is no 
such topological barrier, since the soliton decays to zero at both extrem- 
ities. It can therefore transfer all its energy to the defect mode and cease 
to exist. No soliton-like structure need persist to absorb energy from the 
defect mode. 

This points to a criterion we have not seen mentioned when evaluating 
the effectiveness of collective coordinate ansatzen. If two modes included 
in an ansatz can become 'far from orthogonal' in some parameter regime, 
in that they may be highly correlated or may be used to represent the 
same information, then collective coordinate methods may give mislead- 
ing results for solutions that approach these parameter regimes. 

(2) For intermediate intensities, and specifically r) — 2, the long time behavior 
consists of a captured defect mode and a transmitted soliton. The faster 
the soliton's initial approach, the less energy is captured by the defect 
mode. Comparing Figure 4.4 with the lower graph of Figure 6.3, we see 
this 'capture efficiency' phenomenon in both cases. In neither case does 
a critical velocity separate captured from transmitted solitons. 

(3) For small intensities, the soliton is split into three parts: reflected, trans- 
mitted, and captured components. A prediction is made in [6] about the 
amount of energy in each of these three modes, although no comparison 
is made in numerics. Since the ansatz (5.2-5.3) permits only two modes, it 
cannot possibly capture all of this behavior. A three-mode ansatz, includ- 
ing two solitons and a defect mode, was constructed, but not found to be 
useful in studying this system. Nonetheless at very small (resp. large) ini- 
tial velocities, the soliton is almost entirely reflected (resp. transmitted), 
so the two-mode ansatz is reasonable. In these cases, the initial condition 
hes squarely to one side or the other of the stable manifold illustrated in 
Figure 7.4, and the ODE and PDE simulations compare reasonably well. 
At intermediate velocities, the solution reaches a state shown in Fig- 
ure 8.1, which is a reconstruction of the ansatz solution from the ODE 
parameters. In this flgure, the solution appears as a soliton cleaved in 
two by a defect mode. In the full PDE, the two halves of the sohton com- 
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Fig. 8.1. Reconstruction of the spatial structure of solutions from the two-mode 
ansatz (5.2-5.3), compared with full PDE solutions, for orbits reflected and trans- 
mitted in the case of rj = 0.5. 

ponent would separate and proceed in opposite directions. In the ODE 
reduction, they are unable to do that: the ansatz (5.2-5.3) constrains 
recombination into a single soliton. 

These comments illustrate what we believe to be rather general issues relevant 
to understanding the successes and failures of finite-dimensional or collective 
coordinate ODE reductions in reproducing PDE dynamics. Of course any cor- 
respondence between the solutions to a PDE and its variational ODE approx- 
imation depends on the assumption that the PDE solution remains close to 
the ansatz used in the approximation. This may or may not happen, and it 
is risky to draw quantitative information from the ODE model, for example 
regarding how solutions of the PDE depend on a certain parameter that is 
varied. Indeed, many such studies amount to numerically solving both the 
ODE and PDE, and noticing that the behavior is similar. 

A more nuanced approach that is advocated in this paper and others, is that 
the ODE can be used to illuminate a mechanism that underlies a behavior 
seen in numerical simulations of the PDE. For example, we have seen that 
the existence of a critical velocity in the PDE can be explained by finding 
separatrices in the reduced ODE dynamics. 

We may even take this a step further. An eventual goal of this research is 
to understand the behavior of gap solitons interacting with defects in Bragg 
gratings [2]. The derivation of a variational ODE for that system is compli- 
cated by the fact that gap solitons possess internal modes, which would require 
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additional degrees of freedom in any collective coordinate ansatz, as in [34]. 
Nonetheless, NLCME gap solitons interacting with defects share many qual- 
itative behaviors with the ODEs (5.9) derived in this study. We may draw a 
bifurcation diagram much like Figure 2.1 for this PDE. Further, low ampli- 
tude gap solitons are either coherently transmitted or reflected, whereas high 
amplitude gap solitons are captured when sufficiently slow, and pass by the 
defect if they have enough kinetic energy. Both of these behaviors were seen in 
the numerical experiments of Section 6. We may postulate that mechanisms 
similar to those seen in section 7 are responsible for the "resonance" between 
solitons and defects described in [2] . 
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Appendix: Detailed calculations 

Fixed points on Vq 

We verify the claim that the fixed points (■(/', a) = (vr, a*) on Vq are saddle 
centers. The center behavior within Vq follows from the structure of the re- 
stricted Hamiltonian (7.4), and behavior transverse to Vq is determined by 
the hnearization given in (7.15), evaluated at (tt, a*). The resulting (constant) 
matrix has zero trace and determinant 

D = (c - a*f ({c - a*) - y^^^j^^l) tl^ (c - a*fb{a*), (8.1) 

so, provided D{a*) < 0, the remaining eigenvalues A = ±\/—D are real, 
implying hyperbolic saddle type behavior in the {Z, V) directions. Note that 
D{a*) < iff 

a*>a,^ (8.2) 
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Prom (7.13) the fixed point value a* is given by solution of 

r/ ^ def , / , (2a^ - CO - 1) def x o\ 

^(«) = ^ + (c+l)a= — ^Tp^f — ^ R{a). (8.3) 

Now L(a) is monotonically decreasing and R{a) is monotonically increasing 
(in the range a G (1, c)), so if we can show that L{ac) > R{ac), it follows that 
a* > ttc and hence that D{a*) < 0, as required. But 

(c2 _ c - 1) 1 
L{ac) = 2^ and R{ac) = --, 

and so the claim is true. 



The cross section 



To verify that — {{Z, V,ip,a)\a E {a*,c),ip = tt} is a cross section for the 
fiow it suffices to show that 7^ on E^r. Prom (5.9c) we have 



c^-2ca y2 {2a?-ca-l) , ^ 

V'lv=7r = ^ + (c - a) sech Z , ^ ^ — sech Z 

^ ^ Va — 1 



^2 

= — ^ + (1 - sech Z) 



(? — 2ca 



< (1 - sech Z) 



■(2a2 - CO - 1) 

- 2m 



;2 - 2ca 



c — a) sech Z 
+ (c - a) 



sech Z 



(c — a) sech^ Z 



— [i?(a) — L[a)] sech Z, 



(8.4) 



where R{a) and L(a) are defined in (8.3). Now sechZ < 1 for |Z| < 00, so 
the sign of the leading term in (8.4) is determined by the expression in square 
brackets. The second term of this is always negative for 1 < a < c and the 
first is also negative for a > a*, since > Oc > c/2, as shown above. Finally, 
since R{a) is monotonically increasing and L(a) monotonically decreasing and 
a > a*, the point at which R{a*) — L{a*), the last term is also strictly negative. 
We conclude that ^0 < on E^, as required. 



Reduction and the Melnikov function 



We summmarise the modified reduction procedure and Melnikov calculation 
developed by Holmes and Marsden [35] for two degree-of-freedom Hamiltonian 
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systems in the form (7.7), in which the frequency of the action- angle mode 
depends upon the phase variables in the other degree of freedom. (The pro- 
cedure is also outlined in [36], where it is applied to Kirchhoff's equations for 
equilibria of an elastic rod.) As in Melnikov's 'standard' method [26,27], trans- 
verse intersections of stable and unstable manifolds of a perturbed system are 
found by examining the zeros of an integral computed along the homoclinic 
orbit of the unperturbed system. 

Consider the perturbed two-degree-of-freedom system with Hamiltonian 

H = Ho{V,p; I) + i^Hi{V,p; 1,^1:) = h (= const.), (8.5) 

and let 

r)M f)J-f 

^-^-^ + Oil^) = ^0 + Oi„). (8.6) 

Then, provided ^2 > 0, Equation (8.5) may be inverted and solved for / in 
terms of V, p, ip and the constant h. As shown in [35], we may therefore 
eliminate / and replace time by the conjugate variable and write the reduced 
three-dimensional system on the constant energy surface as a periodically 
forced single degree of freedom system with Hamiltonian —I{V,p,iJj]L,h), 
and evolution equations 

where (•)' denotes d/dip{-). This implies conservation of three-dimensional 

phase-space volume on the constant energy manifolds, and area preservation 
in the two-dimensional Poincare maps defined below. Moreover, in [35] it is 
shown that the reduced Hamiltonian / of (8.7) may be written 

I = Io{V,p; L, h) + /xTi(y,p, V'; L, h) + 0{^x''), (8.8) 

where Xi is 27r-periodic in ^0. Thus, reduction yields the standard form of 
a periodically perturbed single-degree-of freedom system for application of 
Melnikov's method [26,27]. In fact, inserting the series (8.8) into the Hamil- 
tonian (7.7), we find 

XQ^Ho{y,p)-\h), (8.9) 
_ H,{V,pM) 



When fJ, — 0, the reduced Hamiltonian system (8.7) has a phase portrait which 
coincides with that of the full system, since the vector field is given by 

( dl dl\ 1 (dH dH\ 



dp ' dV J Qq \ dp^ dV 
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it therefore also has a homochnic orbit. When > 0, the system (8.7) is non- 
autonomous, and thus we may no longer draw a phase portrait, but we may 
instead construct the Poincare map on the cross section E^^ — {{V,p,ip — 
t/jo)} [27]. By a theorem of McGehee [25] the periodic orbit at infinity, /5, and 
its stable and unstable manifolds W^{P) and W'^{P) persist for small values 

of fjL. 

To prove transverse intersection of W'^{P) and W^{P), we apply the Melnikov 
method to the Poincare map P^g : E^^ —>■ E^^ that results from following the 
fiow from ijj — ijjQ to ip = ipo + 27r. As noted at the beginning of Section 7, and 
treated in greater detail for the analogous Sine- Gordon problem in [5], a result 
of McGehee [25] allows us to apply the Melnikov method even though the fixed 
point at infinity is not hyperbolic. The Melnikov integral can be interpreted as 
a normalized distance between the stable and unstable manifolds at a specified 
point on the cross section E. As in [36], we may then apply a version of the 
usual Melnikov method [26,27] to the reduced system (8.7). This would lead to 
a Poisson bracket involving involving Xq and Xi in the Melnikov integrand, but, 
via Equations (8.9-8.10), this is equivalent to the 'p — V Poisson bracket of the 
original functions Hq and Hi. This yields Theorem 2 as stated in Section 7.1. 
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